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Infrared Behaviour and Running Couplings in Interpolating Gauges in QCD 
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We consider the class of gauges that interpolates between Landau- and Coulomb-gauge QCD, 
and show the non-renormalisation of the two independent ghost-gluon vertices. This implies the 
existence of two RG-invariant running couplings, one of which is interpreted as an RG-invariant 
gauge parameter. We also present the asymptotic infrared limit of solutions of the Dyson-Schwinger 
equations in interpolating gauges. The infrared critical exponents of these solutions as well as the 
resulting infrared fixed point of one of the couplings are independent of the gauge parameter. This 
coupling also has a fixed point in the Coulomb gauge limit and constitutes a second invariant charge 
besides the well known colour-Coulomb potential. 
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I. INTRODUCTION 

The precise mechanism responsible for the confinement 
of the coloured states of QCD is still elusive. Lattice 
calculations of the gauge-invariant Wilson loop have pro- 
vided clear evidence for a linear rising potential between 
static colour charges. The underlying long range inter- 
action, however, is provided by gauge dependent objects. 
It thus seems possible that the confinement mechanism 
itself looks different, depending on the gauge one is study- 
ing. Indeed, recent investigations of the infrared be- 
haviour of QCD in Coulomb and in Landau gauge sup- 
port this picture. 

In Coulomb gauge the confinement of two static exter- 
nal colour charges is related to the colour-Coulomb po- 
tential V cou \(R) 0,0- This quantity is a renormalisation- 
group invariant and couples universally to colour charge. 
It has been demonstrated that the potential is long 
ranged and provides an upper bound for the Wilson po- 
tential V(R) |3- Thus a necessary condition for the Wil- 
son potential to be confining is that V cou \(R) be confining. 

In terms of Green's functions, the colour-Coulomb po- 
tential is given by the instantaneous part of the 00- 
component of the dressed gluon propagator in minimal 
Coulomb gauge, 



g 2 D o(x,x ) = (gA%(x,x o )gA b (0,0)) 

= V cou i(\x\)S(x ) + (non-inst.). 



(1) 



It is a renormalisation group invariant and can be calcu- 
lated explicitly via 3 



V coul (x - y)S ab = fdM-^AK-^M- 



(2) 



where M(A) = —diDi(A) is the Faddeev-Popov oper- 
ator, and the gauge-covariant derivative is defined by 



[Di(A)i 



a. 



gf abc u) t A\. The expectation value 



of the inverse Faddeev-Popov operator is the ghost prop- 
agator, iA(x - y)5 ab = (M(A)~ 1 )% b y , and the long range 
properties of the Coulomb-potential are related to the 
infrared behaviour of A(x — y). Moreover V cou \(R) has 



been calculated both analytically p| and numerically [|| , 
and found to be (almost) linearly rising at large R. A 
running coupling a cou \(k) may be introduced by Q 



Vcoul(fc) 



47ra cou i(fc) 



k 2 



llN c -2N f ' 



(3) 



which for a linearly rising V con i(R) is then proportional 
to 1/fc 2 ; a clear instance of the notion of infrared slavery. 

The situation in Landau gauge, on the other hand, 
appears to be much more intricate. Compared to 
Coulomb-gauge there is no quantity that bears an obvi- 
ous correspondence to the colour-Coulomb potential. A 
renormalisation-group invariant running coupling, how- 
ever, can be defined from either of the primitively diver- 
gent vertices Q- A particularly simple example of such 
a coupling is given by the one obtained from the ghost- 
gluon- vertex, 



a(A 2 ) = |-jfc 8 A 2 (k 2 )D(k 2 ), 



(4) 



which involves the ghost propagator A and the scalar 
part D of the gluon propagator, but not the vertex it- 
self i i ng. This fact can easily be traced back to 
the transversal structure of the gluon propagator in Lan- 
dau gauge and will be discussed in more detail below. 
The infrared behaviour of the running coupling I0J has 
been determined analytically, and an infrared fixed point 
at a(0) w 8.915/iV c has been found [13 • On a quali- 
tative level this fixed point behaviour has been demon- 
strated also for the couplings defined from the three- and 
four-gluon vertices Q • Numerical solutions of the Dyson- 
Schwinger (DS) equations for the ghost and gluon prop- 
agators furthermore show that the fixed point is continu- 
ously connected to the well known perturbative coupling 
in the ultraviolet momentum regime 

The aim of this paper is to explore possible connec- 
tions between the Coulomb and Landau gauge results 
described above. To this end we consider a class of 
gauges that interpolates between Landau and Coulomb 
gauge. In section II we will demonstrate the existence of 
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two renormalisation-group invariant running couplings in 
these gauges. One of these couplings will later be shown 
to be an RG-invariant gauge parameter. The other one, 
aj(k), is an analogue of the Landau gauge expression 
0. We derive the Dyson-Schwinger equations for the 
ghost and gluon propagators in the interpolating gauges 
in section III, and provide their asymptotic infrared solu- 
tions in section IV. We find that there is an instability in 
the solution such that, for each value of the interpolating 
gauge parameter a that appears in the local action, there 
is a class of solutions parametrised by a single parameter 
i]. It turns out that ai(k) has a gauge-independent in- 
frared fixed point a/(0) at a value identical to the one in 
Landau gauge. As shown in section V, such a fixed point 
can be also found in Coulomb gauge, although we are not 
able to determine the exact value in this limit. In section 
VI we show that the degeneracy parameter rj that char- 
acterizes the solution of the DSE in a given gauge may 
be identified with an RG-invariant gauge parameter. We 
summarise and discuss our results further in section VII. 



q-i 



FIG. 1: Ghost-gluon vertex DS equation. 



Thus in the interpolating gauge there are two scalar func- 
tions D tr (k) and D^(k) that characterise the gluon prop- 
agator. They and the ghost propagator A(fc) are func- 
tions of the two scalar variables |k| and k , D tT (k) = 
D tT {\k\,k ) etc. 

A useful property of the dressed ghost-gluon vertex 1 
pabc _ y fl (p,q)f abc , in the (transverse) interpolating 
gauges considered here is the factorisation of both the 
incoming and the outgoing ghost momenta, 



(10) 



II. NON-RENORMALISATION OF 
GHOST-GLUON VERTICES AND RUNNING 
COUPLING 

The interpolating gauges we shall consider are specified 
by the (unrenormalised) gauge condition 



ad A + V • A = 0, 



(5) 



with d 1 ^ — (ado, V), and bare gauge parameter < a < 1. 
The values a = 0, 1 correspond to Coulomb- and Landau- 
gauge respectively. By means of an extended BRST- 
formalism, including gauge and Lorentz transformations, 
this class of gauges has been shown to be renormalis- 
able, and physical observables in these gauges have the 
same expectation values as in covariant gauges . This 
class of gauges has also been studied numerically |l4ll5| . 
The partition function Z = J d[A, c, c] exp[— S-pp] is ex- 
pressed in terms of the local Faddeev-Popov action 



S FP (A,c,c) = S Y m+ I d^xiWr^d^Arf+d'c-DiA^}, 

(6) 

where D(A) is the gauge covariant derivative. The class 
of interpolating gauges that we shall be concerned with 
is obtained in the limit £ — > 0. The gluon propagator 
then satisfies the transversality condition 



k'xD^ = ak Don(k) + hD ifl (k) = 0, 



(7) 



which determines the 3-dimensionally longitudinal com- 
ponents of in terms of Dqo, 



Ao(fc) = D 00 (k) 



D tJ (k) 



+ 



k 2 



k 2 k 2 



Doo(fc). 



(8) 



(9) 



where p'^ — (apo,pi) and the function F\ illJ (p,q) is a 
Feynman integral in every order of perturbation the- 
ory. This factorisation of momenta is familiar in Lan- 
dau gauge, a = 1, and holds in the in terp olating gauges 
considered here for the same reason Indeed, fac- 

torisation of the outgoing ghost momentum p' is imme- 
diate, and occurs for any value of the gauge parameter £. 
Factorisation of the incoming ghost momentum q' can 
be seen easily from its Dyson-Schwinger equation, Fig.^ 
and the transversality of the gluon propagator Z? M „ which 
holds only in the limit £ — > 0, V^D^il — q) = q'^D^(l — q). 
Factorisation of both ghost momenta depresses the ul- 
traviolet divergence of T M (<;,p) by a power, so Z\ is fi- 
nite, and may be normalised to Z\ = 1. This implies 
the non-renormalisation of the two ghost-gluon vertices 
g fabdg^a A b c d and g a f abd d Q c a A b c d . (Here Z x represents 
the renormalisation of the pair of these vertices; in a more 
complete notation, each vertex has its own renormalisa- 
tion constant.) 

An important consequence is the possibility of deriv- 
ing expressions for two invariant running couplings that 
are renormalisation-group invariant, one for each of the 
propagator functions D tr and -Doo- It has been shown 
in Ref. [Tj| that the interpolating gauges are multiplica- 



1 The most general form of this vertex contains both colour anti- 
symmetric and symmetric parts. However the colour symmetric 
part does not contribute to the DSE's for the ghost and gluon 
propagators. Indeed, the propagators are proportional to <5 a (,, so 
in these DSE's the symmetric structure constant from the dressed 
vertex gets contracted with the antisymmetric structure constant 
from the other, perturbative, vertex in the ghost loop of the gluon 
DSE and the ghost DSE, and this contraction vanishes. For this 
reason we omitted the symmetric structure constant from the 
start. 
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tively renormalisable, 



A 

9 



A/? Z\ l2 \ A = A ,fl (c, c) = (c H , c R ) 

g R Z g ; a = 



1/2. 



?l/2 



-l/2~l/2 



Z 3 a R , 



11) 



with 7^ for a 7^ 1. The renormalisation of the 
gauge parameter a preserves the transversality relation 
a R d{)A R Q + diAjii = 0. Gauge invariance furthermore 
enforces the Slavnov- Taylor identity 



1 — Z\ — Z g Z 3 Z^ 2 . 



(12) 



(In the notation of 13] Z 3 = Z c Z e , and c is chosen to 
renormalise contragrediently to Ai, like the source Ki, of 
DiC, and the last equation reads simply Z g Z c = 1.) The 
dependence of the renormalised gauge parameter a R (fi) 
on the renormalisation mass /x satisfies the RG flow equa- 
tion, 



dag 



= 7 an, 



(13) 



where 7 = ^fi 31 "^°^ 3 ^ is a power series in g R with 
finite coefficients. 

The spatial gluon propagator renormalises according 

to 



Dij (k, a, A) = D R<ij (k,a R , fx) Z 3 



(14) 



which fixes the renormalisation of the two gluon scalar 
propagators 

D tr (k,a,A) = D%(k,a R ,[i) Z 3 , 
D 00 (k,a,A) = D RfiQ (k,a R ,fi) Z , 

A(fe,o,A) = A R (k, a R , n) Z 3 , 
Ti(p,q,a) = T R ^(p,q,a R ), 

Fo{p,q,a) = T R . Q (p,q,a R ) Z^ 1/2 zl /2 , (15) 

where we have also written the renormalisation of the 
ghost propagator and the two ghost-gluon vertices. Here 
/1 denotes the renormalisation mass, and A is a cutoff 
scale. 

The Slavnov- Taylor identity l|12|) implies that the mul- 
tiplicative renormalisation constants cancel in the two 
products 

g 2 (A) A 2 (fc,a,A) D tr (k,a,A) = 

g R {n) A 2 R {k, a R , n) D%(k, a R , /x), (16) 

3 2 (A) A 2 (fc,a,A) a 2 (A) D oa (k,a,A) = 

9r(v) A|(fc,a fl ,/z) a R (fi) D Rfi0 (k,a R) fi), (17) 

where it is understood that a R = a R {p) runs with /x, as 
does the suppressed argument g R = g R (n), and corre- 
spondingly for a = a(A) and g — g(A). Since the l.h.s. 
is independent of the renormalisation mass /x, the r.h.s. 
is also. These products define two RG-invariant running 



couplings. However because Lorentz invariance is not ex- 
plicit, these products depend on two scalar variables |k| 
and fco- It is useful to define running couplings that are 
dimensionless and that depend on a single scalar variable. 
This may be done in various ways. For later convenience 
we define an RG-invariant running coupling a/(|k|), and 
an RG-invariant running gauge-type parameter a/(|k|) 
by 



aim 



a 2 (|k|) a,(|k|) 



16 g\ 

3 4tt 



dk 
~2~7 



(A 2 fl ^)(|k|,fc ), 



5 32 Or 9r 



w T 



4tt 



J^(A 2 R D Rfi0 )(\k\,k ) 



(18) 



The second invariant vanishes in the Coulomb-gauge 
limit, a R = 0. However it was shown in 0, El that in 
Coulomb gauge an RG-invariant running coupling can be 
defined by 



Qicoull 



|k 2 | f| D^ (\k\ 



(19) 



where £>^ s g (|k|) is the instantaneous part of D Rt oo(k). 
This is consistent with eqs. (|17|I . (I18|) provided that 



a R A R 
urn — 

a^o a A 



(20) 



and this is true if Z 3 ^ 2 Z Q X ^ 2 Z 3 = 1 holds in the 
Coulomb-gauge limit. In fact it was shown 0, El that 
in Coulomb gauge this identity does hold, in addition to 
the Slavnov- Taylor identities (|12f> . We conclude that also 
in Coulomb gauge there are two RG-invariant couplings, 

^coul 

(|k|), previously known, and a second one, a/(|k|), 
defined above, that to our knowledge has not been con- 
sidered before. 



III. DS EQUATIONS 

It is of interest to solve the Dyson-Schwinger (DS) 
equations in the interpolating gauge, represented by 
Fig. 2, in order to connect the Landau and Coulomb 
gauges which, as discussed above, have very different be- 
haviour. The unrcnormalised DS equation for the ghost 
propagator reads 



A _1 (p) =p',j>n- 



0^ J d*k p>^(k) 
xA(p + k)T u (p + k,p), 



(21) 



where p^ = {ap 0l pi). 

A new element arises when writing the DS equation for 
the gluon in interpolating gauge, as compared to Landau 
gauge. We treat the transversality condition on-shell, so 
projectors onto the transverse subspace appear on the 
right hand side. In interpolating gauges the transverse 
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where 



FIG. 2: DS equations for the gluon and ghost propagators. 



subspace is orthogonal to the vector k'^ — (ako, k), so the 
projector onto the transverse subspace depends upon the 
gauge parameter a which gets renormalised. To avoid 
this complication, we multiply the DS equation for the 
gluon on the left and right by D K \(p), and obtain 

D M = 7§&£WP) / d'k k' x A(k) 



(27T) 4 

xT K {k,p + k) A{p + k) D KV {p) 



(22) 



where the ellipses represent the tree-level term and gluon 
loops. 

If one substitutes in renormalised quantities, using 
(TL2l . these equations read 



d 4k PR.u, D R,^( k ) 



f^f^^o (27r) 4 J "~ '" t ' H ^ tl '^' J 

xA R (p + k)T R ^(p + k,p), (23) 

DrMp) = 0^D R ^ x (p) J d 4 k k' RA A R {k) 

xT R . a {k,p + k) A R (p + k) D Ri<7V (p) + ... , (24) 

where p' — (apo,Pi) and p' R — (a R po,Pi). To obtain 
a finite ghost equation, without divergent integrals and 
without divergent renormalisation constants, we use the 
factorisation of the incoming ghost momentum to write 
Pr,„(p + k >P) = Vr,v\(p + k,p) p' RjX , and we add and 
subtract the part of the loop integral that is quadratic 
in p, 



Pf^CfivPv 



v'tiP^- ^% J d A k PR tli D RtliU (k) 

xA R (k)V R ,„x(k,0)p'„ x . (26) 



The subtracted integral in l|25|l is of higher order in 
p. It is also finite due to the factorisation of exter- 
nal ghost momenta, and only finite renormalised quan- 
tities appear in the ghost DSE, so the remaining term 
PfiCfivPv — cip 2 + c-2p\ which is quadratic in p (and rota- 
tionally invariant) must also be finite. 

We recall the discussion of , and refer to it for more 
details. To avoid Gribov copies, the functional integral 
in A-space gets cut off at the (first) Gribov horizon which 
occurs where the Faddeev-Popov operator M has a van- 
ishing eigenvalue. The Faddeev-Popov determinant van- 
ishes on this boundary, and with it the integrand of the 
functional integral. Thus there is no boundary contri- 
bution to the DS equations, which are consequently un- 
changed by the cut-off at the Gribov horizon. 

Moreover it has been shown in lattice gauge theory Jl8| 
that, as a consequence of the cut-off of the functional 
integral at the Gribov horizon, the ghost propagator A(p) 
is more singular than the free propagator. 2 Thus the 
condition that the functional integral be cut-off at the 
Gribov horizon is imposed by choosing a solution A R (p) 
which is more singular than the free propagator l/(p' -p) 
at p = 0. We therefore require that the term quadratic in 
p on the r.h.s. of l|25|) vanish, p^c^ v p v = 0. This gives the 
finite renormalised ghost DSE, with the horizon condition 
imposed, 



&r(p) 



d 4 k p' R uD R ^ v (k) 



(2tt)< 

x[A R (k) V R ,„ x {k,0) 

- A R {p + k) V R , vX {p + k,p)] p' RtX . (27) 



[As a by-product we note that the vanishing of 
PtiC^Pu — cip 2 + c 2 p 2 = on the l.h.s. of JSHJ) gives 
formulas for the renormalisation constants, 



Z, = 



(27T) 



f 



6 i=i 



»(M) Mfc) 

(28) 



7 7 l/2 



Ng R 

(2tt)4 



d A k D Rfiu (k) V R>vQ {k,0) A R (k) or, 

(29) 



A R 1 (P) = Pp. c p.vPv + 



d 4k Pr^Dr,^^) 



(2tt)4 
x[A R (k) V R<vX (k,0) 

- A R (p + k) V RtU x(p + k,p)]p'„ x , (25) 



2 A simple, hand-waving argument why this should be so is that 
the ghost propagator A = (M~ 1 (A)) blows up on the Gribov 
horizon, and entropy favours configurations near it, which leads 
to an enhancement of the ghost propagator in momentum space 
at p = 0. 
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where an ultraviolet cut-off is understood, and we used of variable 

1/2 — 

Z 3 Z% = 1/Zg, The last equation is of interest because 

1 /2 

in Coulomb gauge it was shown that Z g Z = 1.] 



IV. ASYMPTOTIC SOLUTION IN INFRARED 

In Landau gauge the ghost propagator A(fc) is en- 
hanced at small k whereas the gluon propagator D^ v {k) 
is suppressed at small k. Consequently the ghost loop, 
written explicitly in Q24JI. dominates the DS equation 
of the gluon in the infrared, the remaining terms repre- 
sented by the ellipses, including the tree level term, be- 
ing subdominant [3 • We seek a solution where the ghost 
loop is also dominant in the gluon DSE in the infrared for 
interpolating gauges with gauge parameter near a = 1. 
We suppose that the propagators have infrared asymp- 
totic forms that are homogeneous in p, namely for small 
A we have D(Ap) ss (A~ 2 ) 1+KD D(p), for some anomalous 
dimension kjj, and likewise for A(p). More precisely, the 
infrared asymptotic propagators are defined by 



D(p) = lim A 2+2kd D R (Xp) 

A — ^0 

A(p) = lim A 2+2ka A R {Xp). 



(30) 



They satisfy the infrared asymptotic DS equations, 

Ng% 



dm 



D MX (p) J d 4 k k' R X A(k)k' R ^ 

xA(p + k) D av (j>) , (31) 



A-\p) = 



(2tt)4 

x[A(fc) 



d 4 k p' R ^t)^{k) p' H v 

- A(p + k)} . (32) 



We have dropped the subdominant terms represented by 
the ellipses in the gluon equation (|24J) . We have also trun- 
cated the DS equations by replacing the dressed ghost- 
gluon vertex by the tree-level vertex, which has recently 
been shown to be a good approximation [Til llH l2fT | . It is 
easy to verify, using lfTT|) , fT^ and lfT5|) , that these equa- 
tions are form-invariant under a renormalisation trans- 
formation. 

We now exhibit a change of variables that transforms 
the infrared asymptotic DS equations from the interpo- 
lating gauge to Landau gauge. The gauge parameter a R 
appears in the DSE only in the expressions of the type 
PR,n D R-,nv( k ) = PiDiv+PociRDov. The gauge parameter 
may be absorbed by a rescaling of Dq ^, or of the coordi- 
nate po or both. For this purpose we introduce a factori- 
sation of a R — Or/, where r\ is an arbitrary constant in the 
interval < rj < 1, and 9 = a R \ rf~ . We define the ma- 
trices H = diag(l, 1, 1, 7]) and 8 = diag(l, 1, 1, a_R?y _1 ), 
so Pr uDR,iMu{k) — [pHQD R {k)] u . We make the change 



P M = {Hp)u, = (ilPa,P) 
D^(k) = tT 1 (OD(fe)e)^ 
A(k) = A(k), 



(33) 



In the new variables, the transversality condition reads 
PiiD^ v (p) = 0, and the ghost and gluon DS equations 
become 



A-\p) 



gji J d A kp,D, u {k)p u 
x[A(k)-A(p + k)} , (34) 



^D, x {p) j A(/,) 

xA(p + k) k a D av (p) . (35) 



The change of variable has brought the transversality 
condition, the ghost equation, and the ghost-loop term 
in the gluon equation from the interpolating gauge to the 
Landau gauge. 

These equations have been solved, 0,0], with the 
result, D^ip 2 ) = D L (p 2 ) (5 M „ - p^/p 2 ), 



A(p 2 ) 
g 2 D L (p 2 ) 



b /x 2ka 

(p2)l+KA 

c [i 2KD 



(36) 



where k a = (93 - \/l20T)/98 ~ 0.595353, n D = -2kg, 
with dimcnsionless coefficients b and c. There is an in- 
frared fixed point 

a(p - 0) = (p 2 f t- D L {p 2 ) [A(p 2 )} 2 = °£ 



2tt T(3 - 2k)T(3 + k)T(1 + k) 
~3N~ C r 2 (2- k)T(2k) 
8.915/iV c , 



(37) 



where k = ka- 

In terms of the original variables, this solution reads 



A R (k) -» A(fc) 



/J 



(k 2 + 7 1 2 k 2 y+^ 

2k d 



9 R D u (k) 



X] C jl 



9 2 r d r,oo -* 9RD Q0 (k) = 



(k 2 +r 1 2 k 2 ) 1 + K ° 
a R 2 rf c fi 2KD k 2 
(k 2 + r] 2 k^) 2+KD 



> (38) 



where — ► means asymptotic infrared limit defined in l|30|l . 
For all rj in the range < 77 < 1, the infrared asymptotic 
solution (|38|l is self consistent because counting of powers 
of momenta remains the same as in Landau gauge, so the 
terms in the gluon equation that were neglected remain 
subdominant by power counting. 
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The /^-dependence nicely factors out of the RG- 
invariant products l|16[) and i|17[l 



g\ D ta A 2 = 
a 2 R g 2 R D 0Q A 2 = 



1 



77 cb 2 



(k 2 + 7? 2 fc 2 ) 3 4tt 

k 2 if c b 2 
(k 2 + 77 2 fc 2 ) 4 4tt ' 



(39) 



because kjj + 2ka = 0. The RG- invariant running cou- 
plings, eqs. (JTSJ, become in the asymptotic infrared limit 
a/(0) = lim A ^ a(A|k|), a 7 (0) = lim A ^ a(A|k|), 



a/(0) 



aj(0) a/(0) 



16 /" dfc |k| 5 r? c 6 2 

Y 7 IhT (k 2 + ri 2 kiy 

32 /■ dfc |k| 7 ?7 3 ci 

y 7 yr^r^fcf) 1 4^ 



.(40) 



where we have also rescaled kg — > Afc . This gives 

.915 



«i(0) 
oj(0) = //■ 



cP 

47T 



iV,: 



(41) 



The RG-invariant running coupling a/(|k|) goes to the 
same infrared fixed point, 8.915/iV c , in all interpolating 
gauges as in Landau gauge. On the other hand the RG- 
invariant parameter a/(0) = r\ appears to (and, as we 
shall see in section VI, does in fact) provide an RG- 
invariant characterisation of the gauge. Indeed at 77 = 1 
the solutions (|38|) are Lorentz-invariant, and agree with 
the Landau gauge (apart from normalisation), while at 
7; = they are independent of fco which means instan- 
taneous propagation in space-time, corresponding to the 
Coulomb gauge. The surprise emergence from nowhere of 
the RG-invariant parameter 77 is in fact necessary to char- 
acterise the gauge because the renormalised gauge pa- 
rameter an = an(n) is useless for this purpose, being sub- 
ject to arbitrary renormalisation. Indeed under a finite 

1 /2 1/2 

renormalisation A R ^ — > z 3 Arj and A R o — > z A R ^ 
it changes according to a R — > cirz^ 2 z^ 1 ^ 2 . 

The solution (|38|l has some attractive features. In 
the Coulomb-gauge limit 77 — > 0, the ghost propaga- 
tor becomes independent of fco, as it should. The in- 
frared anomalous dimensions ka and kd of the ghost 
and gluon propagators, and the infrared fixed point a/(0) 
are independent of the gauge parameter 77 in the range 
< 77 < 1, and may have some gauge-invariant meaning. 3 
It is somewhat surprising that the transverse gluon prop- 
agator D R (k) is non-analytic at the unphysical gauge- 
dependent point k 2 + 77 2 fc 2 = 0, whereas it is singular at 
the physical point k 2 + k 2 = in every order of perturba- 
tion theory. However according to the gluon DSE, 
the non-analyticity of the gluon propagator is determined 



by the 77-dependent singularity of the ghost propagator 
which gives the dominant term in the infrared. It seems 
that by imposing the horizon condition, which is non- 
perturbative, the solution is forced into a confined phase 
in which Dp\{k) is non-analytic at an unphysical point. 
Note that D R (k) actually vanishes at this point, because 
1 + Kd < 0, so this non-analyticity does not contradict 
the Nielsen identity which states that poles in physical 
channels are gauge-independent. 



V. COULOMB GAUGE LIMIT 

The infrared anomalous dimensions of D\\ and D R ^o 
are independent of 77, and thus remain constant in the 
Coulomb-gauge limit 77 — * 0. However according to 13811 . 
the asymptotic infrared gluon propagators D tr (k) and 
Doo(k) are proportional respectively to 77 and t; 3 a R 2 
which vanish at 77 = 0, while the second is ill defined 
at au — 0. Because l) tr (fc) and -Doo(fc) represent the 
infrared asymptotic limit of D t R (k) and D Rt oo(k), this 
would indicate that the infrared anomalous dimensions 
change discontinuously at 77 = (while the propagators 
D tr (fc) and Doo(^) at finite momenta remain well-defined 
at 77 = 0). 

Recall that the ghost propagator is given by A(x— y) — 
{{M- r ) xy {A)), where M = -d'^D^ is the Faddeev- 
Popov operator. In the Coulomb-gauge limit, the in- 
verse Faddeev-Popov operator becomes local in time, 
(M~ 1 ) xy (A) -> (M^% y (A)S(x -y ), so the ghost prop- 
agator, when transformed to momentum space becomes 
independent of po, A (p) = A(|p|). The DS equation for 
the gluon is then singular. Indeed in Coulomb gauge the 
ko integration that appears in (|24|) is given by 



dko k' x A R (|k|)A H (|k + p|) k'j 



R.a' 



(42) 



a catastrophically divergent integral. This is the famous 
energy divergence of the Coulomb gauge. It is cancelled 
by a corresponding contribution from the gluon loop 
which we have not evaluated because it is subdominant 
in the infrared for 77 > 0. Thus the asymptotic infrared 
solution we have obtained no longer holds at 77 = 0. 

However we get some information from the ghost 
equation. In Coulomb gauge the time components of 
the gluon propagator do not contribute to the ghost 
DS equation with tree-level vertex, p' R D R tliv {k)p' R v = 



PiDR ti j(k)pj 



(p-k) 2 ] D%(k,k ), by ©. With 



Afl(fc) = Afl(|k|) independent of fco in Coulomb gauge, 
the fco-integral in the ghost DSE, gives the equal- 

time gluon propagator, 



(43) 



D£(k) = (27^ / dfco^(k,fc ) 



3 From ref. |21| it seems also that the anomalous dimensions are 
independent of the gauge parameter £ of linear covariant gauges. 



The ghost equation in Coulomb gauge reduces then to 
a purely spatial equation relating functions of a single 
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scalar variable, 



A^dpl) = 



Ng 2 
(2^)3 
x[A fl (|k|) 



d 3 k [p 2 -(p-k) 2 ] D*(|k|) 
A fl (|p + k|)]. (44) 



In the infrared asymptotic limit, the propagators obey 
power laws, 



Ah(IpI) 

9 2 rDH(\ P \) 



b(p 2 ) 



( p 2)l+7A 

IpI c (m 2 )^ 



(p 2 ) 



1+71 



(45) 



where 7a and jd ar e infrared anomalous dimensions that 
characterise the Coulomb gauge. (The extra power of |p| 
reflects the dimension of D^ t .) Counting powers of mo- 
menta in the ghost DSE (|44H gives the relation between 
the anomalous dimensions, 



7-D = -27A- 



(46) 



The running coupling aj (|k|) defined in 118fl is also an 
RG-invariant in Coulomb gauge where it simplifies to 



a/dpi 



16 Or . 
46 2 c 



A 2 (|p|)^ t (| P | 



3?r 



(47) 



To determine the value of the infrared fixed point b 2 c, one 
would have to also solve the gluon equation in Coulomb 
gauge in the asymptotic infrared. However since ai(0) 
is independent of the gauge parameter r\ in interpolating 
gauges, it may possibly have the same value in Coulomb 
gauge, namely 46 2 c/37r w 8.915/7V C . 

The RG-invariant coupling a/(|p|) approaches an in- 
frared fixed point for |p| — > 0. This is in sharp contrast 
to the behaviour of a cou i(|p|), defined in fj3, which has 
been approximately determined in Refs.[<4 I22L l23l| to di- 
verge as 



acoui(p/A Q cD) ~ y^p, 



(48) 



with (3 m 2. 



VI. GAUGE EQUIVALENCE AND GAUGE 
INSTABILITY 



We considered interpolating gauges that are charac- 
terised by a gauge parameter an (or a) that appears in 
the local action. We found that the space-time character 
of the solutions of the DSE in the infrared limit is inde- 
pendent of an, but depends instead on an RG-invariant 
parameter rj. Thus our solutions are characterized by 
a single parameter, just as gauges are characterized by 
a single parameter, but it's a different parameter. This 



should not come as a surprise. Recall that an satisfies the 
RG flow equation (|13[1 with solution an — a^(/i, c), where 
c is an arbitrary constant of integration. Quantities such 
as g R DnA R that are RG-invariants must be independent 
of fi and thus depend on c rather than an- The gauge 
dependence thus shifts from an to c. In our solution, 
the RG-invariant products g R D R A R and a R g R Dn,ooD R , 
given in 1)39(1 . are found to be independent of an and to 
depend only upon rj. This leads us to make the identifi- 
cation c ~ i], which makes 77 a true (RG-invariant) gauge 
parameter. It is very encouraging that our asymptotic 
infrared solutions display exact features of the renormal- 
isation group even though we are far from the perturba- 
tive regime and used truncated DS equations. 

The possibly novel element is that rj appeared in our 
solution of the DSE as a degeneracy parameter that char- 
acterizes different solutions in a fixed gauge, that is, for a 
fixed value of an- This is true even for the Landau gauge 
an = 1, which is generally thought to be a fixed point 
of the RG group. However our discussion shows that the 
Landau gauge is unique only if one restricts solutions of 
the DSE in Landau gauge to those that are manifestly 
Lorentz covariant. If one drops this restriction, as we 
do, then the DSE in Landau gauge also has a continuum 
of solutions, characterized by the i?G-invariant param- 
eter r). Thus we may say that each gauge is unstable, 
having a continuum of solutions, and moreover different 
interpolating gauges, including the Landau gauge, are 
equivalent to each other. 

If the identification of the degeneracy parameter 77 with 
the RG-invariant gauge parameter c is correct, then it is 
to be expected that the degeneracy of solutions that we 
found (for a fixed gauge) holds not only for the infrared 
asymptotic solutions obtained here, but extends to so- 
lutions of the full DS equations at finite momentum k. 
Namely for each value of the parameter an that appears 
in the local action, there is a class of gauge-equivalent so- 
lutions parametrised by the RG-invariant gauge param- 
eter rj, and moreover the dependence on an is trivial. 
Thus for any value of an, including the Landau-gauge 
value an = 1, there is a one-parameter class of solutions 
of the DSE that continuously connects the value rj = 1, 
where Lorentz invariance is manifest, to a solution with 
77 — >• that approaches the Coulomb gauge, where one 
has a simple confinement scenario. The observation that 
the infrared fixed point of the running coupling and the 
infrared anomalous dimensions of gluon and ghost prop- 
agators are all independent of the two gauge parameters 
77 and an, and the observation that the infrared fixed 
point persists in the Coulomb gauge limit, may, hope- 
fully, bring us closer to a unified picture of confinement 
in Landau and Coulomb gauge. 



VII. SUMMARY AND CONCLUSION 

To summarise, we have found two renormalisation- 
jroup invariant running couplings, olj and a 2 a/, in a 
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class of gauges that interpolate between Landau gauge 
and Coulomb gauge. These gauges are conveniently char- 
acterised by the RG-invariant, r\ = a/(0), that varies 
between zero and one. We have determined the infrared 
behaviour of these couplings from a coupled set of Dyson- 
Schwinger equations for the ghost and gluon propaga- 
tor. In all interpolating gauges, < T) < 1, we found 
an infrared fixed point for the running coupling ai(k) 
at aj(0) = 8.915/A. This fixed point as well as the 
infrared anomalous dimensions of the ghost and gluon 
propagators are independent of the gauge parameter r), 
and coincide with the corres pon ding well known values 
in the Landau gauge 

IS Ellll- Thus they may have 
some gauge-invariant meaning. Indeed, the violation of 
the cluster decomposition principle, a necessary condi- 
tion for confinement, entails that long range correlations 
must be present in Yang-Mills theory. One is tempted 
to speculate that these correlations are triggered by the 
fixed point behaviour of the theory in the infrared. This 
is manifest in Coulomb gauge, 77 = 0. Although we can- 
not determine its exact value in the Coulomb gauge limit, 
we still find a fixed point for ai in the infrared. The other 



coupling, ajai, is replaced by the familiar Coulomb cou- 
pling a cou i, eq. (|19f) . which diverges in the infrared and 
determines the static colour Coulomb potential, eq.©. 
Thus the running coupling, that in Landau gauge merely 
displays a fixed point in the infrared, bifurcates by gauge 
instability into two RG-invariant running couplings in 
interpolating gauges, and these in Coulomb gauge are 
responsible for fixed point behaviour on one hand and 
on the other for the long-range colour-Coulomb poten- 
tial between static colour charges. 
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